Scanning tunneling microscopy of adsorbed molecules on metalic surfaces for nearly 

localized atomic states 
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We consider a Hubbard-Anderson model which describes localized orbitals in five different sites 
hybridized both among themselves and with a continuum of extended states. A square planar geom- 
etry with an atom at the center is used to represent TBrPP-Co molecules. When the renormalized 
effective hopping between sites is small compared with a Kondo energy scale determined by the site- 
continuum hybridization, the system can be described as a set of independent Kondo resonances, 
rather than molecular states. We study the crossover between both regimes and analyze the spectral 
density of conduction electrons as a function of position. The results are in qualitative agreement 
with measurements of the differential conductance in a system with TBrPP-Co molecules adsorbed 
on a Cu(lll) surface. 

PACS numbers: 73.22.-f, 75.20.Hr, 68.37.Ef 
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I. INTRODUCTION 

The transport of electrons through molecules is rele- 
vant in many branches of science and important for po- 
tential applications in electronic devices 1 - - — and molecular 
spintronics^* 5 -. On the other hand, the study of many- 
body phenomena in nanoscale systems has attracted 
much attention in recent years. In particular the Kondo 
effect, which arises when a magnetic impurity interacts 
with a continuum of extended states, has been observed 
in a variety of nanoscopic systems 4 ^— . Scanning tun- 
neling spectroscopy (STS) has made it possible to probe 
the local density of states and Fano antiresonances have 
been observed for several magnetic systems on metal 
surfaces 4-9 . These antiresonances observed in the differ- 
ential conductance, reflect a dip in the spectral density 
of conduction states near the Fermi level caused by the 
Kondo effect 11 . The half width of the dip is given by a 
characteristic energy scale, the Kondo temperature Tk- 
Furthermore, corrals built on the (111) surface of noble 
metals or Cu have been used to project the spectral fea- 
tures of the Fano-Kondo antiresonance (FKA) to remote 
place o 7 ' 11 . The observed Fano line shapes for one mag- 
netic impurity on these surfaces have been reproduced 
by many-body calculations 1 ^—. For corrals, predictions 
of the variation of Tk with the position of the impurity 
were mad o 11 ' 18 . 

When several magnetic impurities interacting between 
them are present, the theoretical analysis is more diffi- 
cult. For Cr trimers on the Au(lll) surface, controversial 
results were obtained^— which were in contradiction 
with at least part of the puzzling observed dependence of 
Tk with geometry 9 -. More recently, a qualitative expla- 
nation was provided using a Hubbard- Anderson model 
(HAM) described in Section II2 3 .. 

Recent STS results for TBrPP-Co [5, 10, 15, 20- 
Tetrakis-(4-bromophenyl)-porphyrin-cobalt] molecules 



adsorbed on a Cu(lll) surface are also surprising 5 -. The 
width of the FKA depends on the position of the mi- 
croscope tip along the molecule, varying within a factor 
two among different positions 24 . Quite generally, for 
any system in which only one atom is highly correlated 
(say Co), if the spectral density of the system without 
this magnetic impurity is featureless around the Fermi 
energy ep, the STS near ep is dominated by the Green's 
function of the impurity and the width of the FKA is the 
same everywhere. Only the intensity varies with position 
(Eqs. (11) to (14) of Reffii). Experimentally, for a 
sysyem with one impurity at the focus of an elliptical 
quantum corral, the observed FKA is similar at both 
foci and only the amplitude differs 7 In systems with 
several magnetic atoms with a sizable hopping between 
them, one can define an effective Hamiltonian in terms 
of molecular extended orbitals and if Kondo physics is 
present, again the energy dependence of STS near ep is 
essentially independent of positione r 23 ' 25 . This results 
are confirmed by our calculations on the HAM. 

In Ref£, the experimental results were interpreted in 
terms of a single impurity Kondo model, with a position 
dependent exchange interaction which increases linearly 
with the local spin density. This is strictly valid only in 
the limit in which the system can be considered as a con- 
tinuum of independent magnetic impurities. This situa- 
tion resembles the Kondo lattice in the limit in which Tk 
is larger than the effective intersite exchange J 26 . How- 
ever, in the molecular system, the competing parameter 
is the magnitude of the intersite hoppings t (rather than 
7), neglected in the simplified analysis of RefA Thus, 
one expects a crossover between a regime with indepen- 
dent "impurities" for small t to a collective behavior for 
large t. In the limit of an infinite molecule disconnected 
to the metal, one expects that this crossover turns to a 
Mott localization transition 27 . 

In this work, we describe the system with a HAM in 
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which each magnetic site is hybridized with a continuum 
of extended states and in addition there is a hopping 
between magnetic sites. We obtain that for reasonable 
parameters (most of them used before) , the essential fea- 
tures of the experiment of Ref£ are reproduced. More- 
over we study the above mentioned crossover between 
localized and extended behavior. 

In Section II, we describe the model. The STS spectra, 
comparison with experiment, and the crossover from the 
regime of nearly isolated impurities to clllective molecu- 
lar behavior (either with one Koddo peak or with split 
Kondo peaks) is presented in Section III. Section IV con- 
tains a short summary and discussion. 



II. MODEL, APPROXIMATIONS AND 
PARAMETERS 

The Hamiltonian can be written aa 2 ^ 

H = H mo \ + H mct + Hy, 
H mo \ = ^ Eid\ a d la + Uidl^di^dl^dii 

- ty(4r«**r+H.c.) 1 

Hmet — ^ ekC ko .Ck(rj 

H v = J2(y^ ikRs dlc ka + R.c). (1) 
kjtr 

H mo \ is a Hubbard model that describes the isolated 
molecule. Here d\ a creates an electron at site i lying 
at Ri with spin a. H mct corresponds to the conduction 
band of the metallic substrate, and Hy is the hybridiza- 
tion between the atoms of the molecule and the extended 
states of the metal, assuming that the wave functions of 
the latter are plane wavesii. Following earlier work 28,29 
we represent the molecule with 5 sites, one at the center 
representing the Co atom and four at the corners of a 
square representing bromophcnyl lobes. The model has 
Ci v symmetry. We assume that the tip of the scanning 
tunneling microscope (STM) is far enough from the sys- 
tem in such a way that it does not alter its electronic 
structure. 

This model is similar to the one proposed by Chiappe 
and Verges for similar molecules 2 ^* 2 ^. The main differ- 
ence is that we consider that not only Co, but all sites 
are highly correlated, as suggested by recent ab initio 
calculations 5 . In addition Hy is somewhat more realistic 
and contains phase factors exp(zk ■ Rj), which were es- 
sential in the description of Cr trimers on Au(lll)^. We 
also take Ui — > +oo. Since the Ui are expected to be of 
several eV, the Uj were proposed of the order of 0.2 eV, 
and the resonant level widths T^pVf are smaller than 1 eV 
(see below), this assumption is not a serious one, and on 
the other hand allows us to use the slave-boson mean-field 
approximation (SBMFA) introduced by Colema n 25 i 30 ' 31 . 



In this treatment, the creation operators for local- 
ized electrons are represented as a product of a boson 
and a new fermion as d\ a — bifj^, with the constraint 
b\bi + J2a flafia = !• The approximation consists in re- 
placing the boson operators b\, bi are replaced by their 
expectation values (b\) = (bi)* obtained minimizing the 
free energy. This approximation reproduces correctly the 
spectral density near the Fermi level, and the exponen- 
tial dependence of Tk with the parameters in the Kondo 
regime. 

As in a previous work^ 3 -, we have approximated the 
angular average of J2k ex P[ik ' (R-i ~~ Rj)] > by its value at 
the Fermi surface fyj . For an isotropic three-dimensional 
(3D) band /3ij = sxa(kprij)/(kprij), where k,F is the 
Fermi wave vector and r.y = |Rj — Rj|. We take 
k F R = 3.8, where R is the distance between the Co 
site at the center and any lobe at the corners^ 3 . We 
also assume for simplicity a constant density of conduc- 
tion states p = 0.1/eV. The results practically do not 
change for other values of p if the hybridizations Vj are 
scaled keeping pV^ constant. Following Ref. — , we take 
to = 0.15 eV for the Co-lobe hopping and t\ = 0.2 eV for 
the lobe-lobe one. From experiment^, the first occupied 
orbital is located 0.7 eV below the Fermi level ep and is 
assigned to the 3d z 2 orbital of the Co atom. Choosing 
the origin of energies at ep = 0, then Eq = — 0.7 eV and 
the on-site energy for the lobes should be smaller. We 
take Ei = —3 eV. The hybridizations Vbo = 0.68 eV and 
Vi = 1.68 eV were chosen to reproduce approximately 
the observed width of the resonances when the tip is on 
a specific site. 



III. RESULTS FOR THE TUNNELING 
SPECTRA 

The resulting values of the expectation values of the 
boson operators (&Co) = 0.18 and (6/) = 0.13 indicate a 
strong renormalization of the quasiparticles at the Fermi 
energy. In particular, the effective to is renormalized to 
to{bco)(bco) — 3.6 meV. This is smaller than the re- 
ported Tk- This means that the system is in fact in a 
regime of rather independent "impurities" with two dif- 
ferent Tk- While this picture is useful to gain physical 
insight, the intersite hopping and the iteractions between 
sites through the conduction band are still important and 
affect the width and shape of the STS features, as ex- 
plained below. 

The fact that the renormalization is stronger for the 
lobes ((&;) < (bco)) is consistent with results of ab ini- 
tio calculations which find a larger spin density at the 
bromophenyl units . 

For weak coupling of the STM tip to the system, the 
measured differential conductance dl/dV is proportional 
to the density of states of a mixed operator h a which 
depends on the position R of the tip. For a voltage dif- 
ference V between the metal and the tip one ha s 11 ' 17 
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FIG. 1. (Color online) Differential conductance as a function 
of voltage when the tip is on Co (dashed line) or on a lobe 
(full line) for q Co = 0.136, qi = 0.114. 



FIG. 2. (Color online) Differential conductance as a function 
of voltage for qj = and different positions of the tip along 
the Co-lobe path. 



dl/dV ~ Ph {eV), 

Ph{u) = — (((M^)Ww - {{h a ;h\)) u -ie) , 

MR) = elk ' Rjc k- + &( R K>- (2) 

k j 

where qj (R) is proportional to the hopping between the 
tip and the site j and controls the line shape and its 
asymmetr y 1 1 > 1 7 . It is only significant when the tip is 
right above the site j. Then, here we assume qj(R) = 
qj<5(R" — Rj"), where the superscript on the position 
vectors denote the component parallel to the Cu(lll) 
surface. 

The resulting dl / dV when the tip is either on Co or on 
the lobes is shown in Fig. [U The values of qj were chosen 
to result in a line shape that agrees with experiment^. In 
any case, these values are small indicating that dl/dV 
is dominated by the conduction density of states, which 
presents a Fano-Kondo antiresonance, as for Co impuri- 
ties on Cu or noble metal surfaces^. A fit of the Kondo 
temperature using the simple expression given in ReL s 
leads to T K = 112 K on Co and T K = 208 K on the 
lobes. The main features of the experiment are repro- 
duced. When the tip is on the lobes, we obtain a struc- 
ture near -10 meV. While the experimental resuls seem 
rather noisy, a similar structure seems to be present only 
in some of the observed spectra (for example T— 109 K, 
185 K, 212 K in Fig. 3 (b) of Ref. 5). In any case, 
this structure disappears if a broadening of the quasi- 
particles at finite energies (not taken into account in the 
SBMFA)22, or slightly smaller tij or larger Vi are used. 

In Fig. [2] we show the local conduction density of 
states ph(eV) for qj = (proportional to dl/dV) as a 
function of position along a path between Co and a lobe 
R = aRi + (1 — a)R Co . Starting at the Co site (a = 0), 
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FIG. 3. (Color online) Differential conductance as a function 
of voltage for qj = and different Co-lobe hopping to with 
the tip on top of the Co atom. 



the dip broadens and loses intensity. The intensity is 
recovered at the position of the lobe while the peak con- 
tinues to broaden. The space dependence of the depth 
of the structure might be expected from what we know 
from the one-impurity case, where the amplitude of the 
dip decays as 1/d 2 (1/d) with the distance to the impu- 
rity d for a 3D (2D) conduction band£. Therefore, the 
amplitude of the dip has a smaller variation if a 2D band 
or smaller kp is assumed. In particular, for kp = 0, the 
dependence of the conduction density of states near ep 
with position disappears. 

To study the crossover to the regime in which molec- 
ular states dominate the physics, we have increased the 
Co-lobe hopping to. The result for dl/dV when the tip 
is on top of the Co atom is shown in Fig. [3] For t = a 
fit of the curve leads to Tk = 33 K. This agrees with the 
fact that for an isolated Co atom on a Cu(lll) surface, 
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FIG. 4. (Color online) Differential conductance as a function 
of voltage for qj = and different Fermi wave vectors with 
the tip on top of the Co atom. 



the observed Tk (Ref. 8 ) is about a factor 1/2 of that 
in the molecule (in our case however, there is a remain- 
ing interaction with the lobes through the bulk states 
unless )cf — t oo). As to increases, the peak near the 
Fermi energy splits in two for large enough to ~ 0-2 eV, 
and the splitting increases with t . The physics of this 
splitting is also present in previous experimenta l 33 ' 34 and 
theoretica l 31 ' 35 " — studies of two quantum dots, and for 
impurities with hopping to discrete extended states (as 
surface states inside a corral with hard wall a 13 ' 14 ). This 
structure is also apparent when the tip is on the lobes 
(not shown), but an additional structure at ep remains 
due to states that do not mix with the Co ones for sym- 
metry reasons. This structure tends to split when the 
molecule is distorted to a rectangular shape. However, 
the effect of this distortion is not dramatic if the other 
parameters are kept constant. Nevertheless, it is known 
that this distortion alters the Co-surface distance and as 
a consequence Vq should be modified 4 -. 

For large to and positive Ei , we have verified that there 
is a Kondo effect with a molecular state, and dl/dV has a 



similar shape near the Fermi level for any position of the 
tip, with the same Tk within less than 1%. This situation 
is analogous to the physics of quantum corral o 7 ' 11 ' 13 ' 14 
and Cr trimers on Au(lll)2ii2r— . 

The value of kp affects the interference effects medi- 
ated by the conduction band, between different sites of 
the molecule. For kp — > 0, all conduction states near the 
Fermi energy mix with the same molecular orbital regard- 
less of its direction resulting in a maximum interference, 
For k p — y oo the interference occurs only through the 
hoppings tj. In Fig. [4] we show how the conduction elec- 
tron density of states at the Co site changes with kp. For 
the smaller values investigated, there is a splitting of the 
Fano-Kondo antiresonance due to an effective mixing of 
the molecular states mediated by the band. 

IV. SUMMARY AND DISCUSSION 

For a molecule with strongly correlated states near the 
Fermi energy on top of a metallic surface, we have stud- 
ied the competition between interatomic hopping and 
the hopping of these atoms with the continuum of ex- 
tended metallic states. We show that strong correlations 
can lead to the rather unexpected situation in which the 
hopping to the metallic substrate dominates, leading to 
a physics characterized by different Kondo impurities as 
a first approximation. The resulting dl/dV and its de- 
pendence with position agrees with recent measurements 
on TBrPP-Co molecules adsorbed on a Cu(lll) surface. 
We have also analyzed the crossover to the most usual 
regime in which molecular extended states dominate the 
physics. Our approximation (slave bosons in the mean- 
field approximation) is known to describe qualitatively 
the physics near the Fermi level in both limit s 11 ' 25 ' 31 ' 38 
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